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Abstract 

r~| I In relativistic mechanics the energy-momentum of a free point mass moving without 

acceleration forms a four-vector. Einstein's celebrated energy-mass relation E = m(? 
is commonly derived from that fact. By contrast, in Newtonian mechanics the mass 
is introduced for an accelerated motion as a measure of inertia. In this paper we 
rigorously derive the relativistic point mechanics and Einstein's energy-mass relation 
using our recently introduced neoclassical field theory where a charge is not a point 
but a distribution. We show that both the approaches to the definition of mass are 
^ \ complementary within the framework of our field theory. This theory also predicts a 

(— I ■ small difference between the electron rest mass relevant to the Penning trap experiments 

CLf and its mass relevant to spectroscopic measurements. 

1 Introduction 

The concept of a point particle is at the very heart of Newtonian mechanics. Regarding 
this, Einstein said, "Physical events, in Newton's view, are to be regarded as the motions, 
governed by fixed laws, of material points in space. The material point is our only mode 
of representing reality when dealing with changes taking place in it" |EinIOj . Even in the 
quantum mechanics the concept of a point particle continues to be fundamentally exact. 
As Feynman puts it, "The wave function if) (r) for an electron in an atom does not, then, 
describe a smeared-out electron with a smooth charge density. The electron is either here, 
or there, or somewhere else, but wherever it is, it is a point charge" |FeyIII| . One of our key 
motivations for introducing a neoclassical field theory of distributed elementary charges in 
|BF4] - jBF7] was a desire to account for particle properties as well as for wave phenomena in 
a single mathematically sound Lagrangian relativistic theory. This theory is self-contained 
and, consequently, all particle properties must naturally come out from the field equations as 
approximations. We have already demonstrated that the theory implies in the non-relativistic 
limit (i) the non-relativistic particle mechanics governed by the Newton equations with the 
Lorentz forces and (ii) the frequency spectrum for hydrogenic atoms. In this paper we study 
relativistic aspects of the theory, namely, the motion of a localized wave which describes the 
distributed charge in an external EM field. An idea of a particle emerging as a well localized 
field was explored by a number of authors, in particular in the form of "extended charge" 
models. The Lorentz-Abraham model and its developments was studied and advanced in 
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|AppKie| , jBamGalj . [IKM] . |Ji^ . [Kii2] . |KKSpo|, [M], jPeal], |Roh], |Schwin], [Spo] , 



|Yag| ■ In 1905-1906 Poincare suggested |Poi] (see also [Jack] . |Rohj . |PauRT] . |Schwin] . |Yag 

and references within) adding non-electromagnetic cohesive forces to the Lorentz-Abraham 
model. Here we study dynamics in an external field of a distributed charge with such cohesive 
forces. 

Recall now the fundamentals of the relativistic dynamics of a mass point and the rela- 
tivistic field theory. The relativistic dynamics of an accelerating mass point charge in the 
case where the acceleration is caused by the electromagnetic (EM) external field Eex, Bex is 
described (see, for example, |Barut j . |PauRTj ) by the following equation: 

- (Mv) = fLor (t, r) , M = mo7, 7 = f 1 - ^ 

where v = dr/dt is its velocity, nio is the rest mass of the mass point, fLor is the Lorentz 
force, 

fLor (t, r) = q (Eex (t, r) + -V X Bex (t, r) ) , (2) 



q is its charge and 7 is known as the Lorentz factor. Equations ([T]) for the space components 
of the 4-vector are usually complemented with the time component 

^(MC^) =fLor-V. (3) 

In particular, for small velocities when |v|/c ^ 1 we readily recover from equation ([1]) 
as its non-relativistic approximation Newton's equation with the Lorentz force by setting 
7 = 1. Note that in ([T]) the rest mass mo of a point is an intrinsic property of a point and is 
prescribed. 

In a relativistic field theory the relativistic field dynamics is derived from a relativistic 
covariant Lagrangian. The field equations, total energy, momentum, forces and their densities 
are naturally defined in terms of the Lagrangian for both closed and non closed (with external 
forces) systems (see, for instance, [And] . [Barutj . [Molj . |PauRTj . |Sexl] . see also Section ITTi) . 
For a closed system the total momentum has a simple form P = Mv with a constant velocity 
V, and the relativistic invariance of the energy-momentum 4-vector allows to derive Einstein's 
celebrated energy-mass relation, 

E = Mc^ M = mo7, (4) 

between the mass M and the energy E (see, for instance, |And| Sec. 7.1-7.5], |Moll Sec. 3.1- 
3.3, 3.5], IPauRTl Sec. 37], jSe^ Sec. 4.1]). As stated by Pauh, [PauRTl p. 123]: "We can 
thus consider it as proved that the relativity principle, in conjunction with the momentum 
and energy conservation laws, leads to the fundamental principle of the equivalence of mass 
and (any kind of) energy. We may consider this principle (as was done by Einstein) as the 
most important of the results of the theory of special relativity." Yet another important 
point about the Einstein mass-energy relations is made by Laue [Schilpp , p. 529]: "... we 



can determine the total amount of energy in a body from its mass. We thereby get rid of the 
arbitrariness of the zero point of energy which the former definition of energy (cf. Section 
III) was forced to introduce. There are not merely energy differences, as before; the energy 
possesses a physically meaningful absolute value." 
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Now let us take a brief look at the difference in determining the mass in the dynamics of 
a single point versus a general Lagrangian Lorentz invariant field theory. 

First, since equation ([1]) for a single mass point has a form of Newton's law, one can 
determine as in Newtonian mechanics the mass M as a measure of inertia from the known 
force fLor and acceleration d^r (the variability of 7 can be ignored for mild accelerations). 
In a relativistic field theory for a closed system the energy-momentum is a four-vector, and 
that allows to define the total mass and the rest mass of the system in terms of the energy 
by Einstein's formula dl]) in the case of uniform motion. However, in the case of a general 
non-closed system (which is the subject of our primary interest since we study field regimes 
with acceleration), there is no canonical way to determine the mass, position, velocity, and 
acceleration. For a non-closed system there is even a problem with a sound definition of the 
center of mass since "the centre of mass loses its physical importance" |Molt p. 203]. To 
summarize, in a general relativistic field theory the rest mass is defined for a uniform motion, 
whereas in the Newtonian mechanics the concept of inertial mass is introduced through an 
accelerated motion. 

Our manifestly relativistic Lagrangian field theory describes a single charge by a complex- 
valued scalar field (wave function) ip [t, x) satisfying the nonlinear Klein-Gordon (KG) equa- 
tion, 

- + VV - G' (1^1^) ^ - = 0, (5) 

c X 

where m is a positive mass parameter, and x is a constant which coincides with (or is close 
to) the Planck constant h. The expressions for the covariant derivatives in ([5]) are 

dt = dt + \,^, V = V-^Ae., (6) 
X Xc 

where q is the value of the charge, and (p^^, Aex are the potentials of the external EM field 
(which can be thought of as the one produced by all remaining charges of the original system of 
many charges, see [BF4j - [BF6] for details). The nonlinear term G' is given by the logarithmic 
expression 

G' (s) = (s) = -a-^ [In (ah) + In n^/^ + 3] , s > 0, (7) 

where a is the charge size parameter. It is established below that for certain regimes of 
accelerated motion the relativistic mass point equations ([I])-(l3]) are an approximation which 
describes the behavior of the field ip when it is well localized. The charge localization is 
facilitated by the Poincare-like cohesive forces associated with the nonlinearity G'^{s). In 
particular, a free resting charge with the minimal energy and size a has a Gaussian shape, 
namely, \ip\ = 7r~'^/^a~^/^e~'^' " The localization can be described by the ratio a/Rf 
where i?/ is a typical length scale of the spatial variation of the external EM forces. 

As it is commonly done, we start with assigning a position to the distributed charge by 
using its energy density u (t,x) to define its total energy E and energy center (or ergocenter) 
r (t) by the formulas 

E (t) = j u (t, x) dx, r{t) = -^j^u (t, x) dx. (8) 

If a field ip satisfies the field equation (|5]), its ergocenter r and energy E satisfy equations 
which can be derived from the conservation laws for the KG equation (jS]). We prove under 
the assumption of localization in the asymptotic limit a/Rf ^ that these equations turn 
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into the relativistic mass point equations ([I]), ([3]). Remarkably, the value of the inertial mass 
determined from the equations by the Newtonian approach coincides exactly with the mass 
given by Einstein's formula (jl]). Of course, a convincing argument for the equivalence of the 
inertial mass and the energy based on the analysis of the charge momentum when it interacts 
with the electromagnetic field has been made by Einstein |Ein05aj . but here the same is 
obtained through a thorough mathematical analysis of a concrete Lagrangian model. 

Since our Lagrangian theory is self-contained, Einstein's energy-mass formula or any 
asymptotic law of motion must be derived within the framework of the theory. Importantly, 
since the limit mass point equations are derived, and not postulated, the resulting rest mass 
is shown to be an integral of motion rather than a prescribed constant. Consequently, the 
rest mass may take different values depending on the state of the field. In particular, in 
addition to the primary Gaussian ground state there is a sequence of rest states with higher 
rest energies and rest masses. The possibility of different rest masses comes from the fact 
that in our theory an elementary charge is not a point but is a distribution described by a 
wave function ip. The charge although elementary has infinitely many degrees of freedom, 
with internal interactions of not electromagnetic origin, contributing to its internal energy. 

Our theory yields a simple expression for the rest mass mo of a charge which differs 
slightly from the mass parameter m, namely, 

mar 

mo = m + --|, 9 

where qq is the reduced Compton wavelength. Evidently, the difference mg — m vanishes 
as ac/a — )■ 0. In the nonrelativistic case treated in |BF4] - jBF6] the nonlinear Schrodinger 
equation follows from the KG equation in the limit c — )■ cxo and a — )■ where a = q^x~^^~^ 
is the Sommerfeld fine structure constant. In this limit the mass parameter m coincides with 
the inertial mass mo in evident agreement with the limit of relation (|9]) with a fixed Bohr 
radius a-Q = ac/a. Observe that our theory provides for complementary interpretations of 
the relativistic and non-relativistic masses which are sometimes considered (see |EriVoy| ) to 
be "rival and contradictory." 

The focus of this paper is on the charge motion in an external EM field which has a mild 
spatial variation. The external EM fields with stronger spatial variation are covered by the 
theory as well but with the use of other techniques. Namely, in the case of the hydrogen 
atom |BF6j . |BF7j . for the electron in the Coulomb field we derive the well known energy 
spectrum involving the Rydberg constant i?oo with the mass parameter m in the place of the 
electron mass. The difference between the mass m entering the spectroscopic data and the 
inertial mass mo as it appears in the Penning trap experiments is discussed in Section 12.51 

In the following sections we analyze relativistic features of the charge emerging from the 
underlying field dynamics. In particular, in Section [2] we derive the mass point equations (II]), 
(|3]) under an assumption that the charge wave function remains localized. In Section |3] we 
show that the localization assumption is consistent with the KG equation (|5]). In particular, 
we describe there a class of external EM fields in which the charge maintains its localization 
in the strongest possible form when it accelerates according to the KG field equation (E]). 
The relativistic theory of many interacting charges is left for another paper. 
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2 Relativistic distributed charge as a particle 



In |BF4j - [BF7] we developed a neoclassical theory for many interacting charges based on 
a relativistic and gauge invariant Lagrangian. We have demonstrated there that this La- 
grangian theory describes electromagnetic interaction in all spatial scales: it accounts for 
at least some quantum phenomena at atomic scale including the frequency spectrum of the 
hydrogen atom, and it accounts for the classical motion of non-relativistic charges when they 
are well separated and locahzed. In this paper we study the relativistic aspects of our theory 
in the case of a single charge. Our general Lagrangian in this special case turns into the 
following relativistic and gauge invariant expression 



^1 w = ^\ -Md:r - ■ V V* - Kiri^ - g (^*^) |> , (lo) 



2m \ c^ 

where ip {t, x) is a complex valued wave function over the space-time continuum and ip* is its 
complex conjugate. In the expression (fTOj) c is the speed of light, 

Ko = . (11) 

X 

The covariant derivatives in ffTOj) are defined by ([6]) and the nonlinearity G [s) is defined by 
the formula 

G{s) =Ga{s) = -a-~^s[\n {ah) +\n7r^/^ + 2] , s > 0, (12) 

where a > is the size parameter. The field equation corresponding to the Lagrangian (ITOl) 
is the nonlinear Klein-Gordon (KG) equation (j5]). 

It is proven in pF4j - [BF6j that in the non-relativistic case the trajectory of the ergocenter 
converges to a corresponding solution of Newton's equation with the Lorentz force as a/ Rf — ?■ 
where i?/ is a typical variation scale of the external EM field. It becomes our goal now 
to derive from the KG field equation (O ) the relativistic law of motion for the ergocenter 
r(t). Namely, we show below that the ergocenter trajectory r (t) converges to a solution of 
the relativistic equation ([1]) as a/i? —)■ if the charge is localized. The derivation is based 
entirely on the analysis of the KG equation (j5]) and corresponding conservation laws. 

2.1 Field symmetric energy-momentum tensor and conservation 
laws 

There are a few popular conventions in setting up coordinates and the metric in the Minkowski 
four dimensional space-time. We pick the one which seems to be dominant nowadays as in 
|Barutt Section 1], |LanLif F| Sections 1.1-1.4, 2], [ItzZubt Section 1-1-1]. The time-space 
four-vector in its contravariant and covariant x^ forms is represented as follows 

x = x>' = = (ct,x) , = 0,1,2,3; (13) 

(14) 











d 
^x^^ 


(ia.,v; 


) ' ^ dXf, 
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with the common convention on the summation of the same indices. The metric tensor gf^,y 
g^^ is defined by 



1 














-1 














-1 














-1 



We also use a common convention for the space 3-vector 

^) =x, 2 = 1,2,3, 



(16) 



X 



/"I 2 'HN 



(17) 



emphasizing notationally by the Latin superscript its difference from 4-vector x^ with the 
Greek superscript. 

For vector field potentials (v^. A) in the KG equation (jS]), we use standard relativistic 
notations with the four- vector potential A^, four- vector current density , and the electro- 
magnetic field F^^: 



A^^ = (^,,, AeJ , = {cQ, J) , F^^- = d^A'^ - d-A^^, 



so that 



where, as always. 



ex 



E = -V<^,, - -dtA, 






-El 


-E2 


-Es 


El 





-B3 


B2 


E2 


B3 





-Bi 


Es 


-B2 


Bi 






B = V X Ao 



:i8) 



(19) 



(20) 



The Lorentz force density of the external EM field F^^ acting on a 4-current is of the 
form 

/'^ = -F'^V^= ( -J-E,pE+-JxB ) . (21) 
c \c c / 

We turn now to the analysis of the Lagrangian (fTOj) . As a consequence of its gauge 
invariance, we obtain the conserved four-current J'^ = (cp, J) where the charge density and 
the current are defined by the following expressions, [BF7j . 



Xq ^ dtiJ 

^ Im — + -ip^^ 

mc^ ip mc^ 

^Im — - — A 



(22) 



m 



mc 



They satisfy the continuity equation 



dtp + V -3 = 



implying the charge conservation 



/ p (t, x) dx = p = q, 
7r3 



(23) 



(24) 



where we choose the constant p to be exactly the charge q as it arises in Coulomb's law, see 
[BF4]-[BFf] for details. 
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There are well understood approaches, due to Belinfante and Rosenfeld, |LanLif F\ Sec. 
32], |Barutt Sec. III. 4], |LanVPM| Sec. 22], to constructing the symmetric energy-momentum 
tensor (EnMT) T^^" from the Lagrangian C based on its invariance with respect to the 
Poincare group. The EnMT T^'^ corresponding to our Lagrangian Ci is constructed by the 
same method as in |BF4j and jBF5j yielding as it does there 



where the covariant derivatives ip'^ are defined by 



ip'" 



8'' + 



c 



dt = dt + 



We 

X 



V = V 



igAe 



(25) 



(26) 



and ip*'^ is the complex conjugate to ip''^. 

We proceed with the interpretation of entries the symmetric EnMT T^'^, |LanLif Fl Sec- 
tion 32], [MorFeshll Chapter 3.4] 



U 



cp 
-a' 



-a 



cp^ 
—a 



-a 



12 
22 



cp^ 
—a 



-a 



13 
23 



I53 _^31 _^32 _^33 



u 


energy density. 


p) 


momentum density. 


= c^p^ 


energy fiux density. 


0'^ = a^' 


symmetric stress tensor. 



(27) 



where i,j = 1, 2, 3. 

Remark 1 As a consequence of the symmetry of EnMT T^'^ we have the following relation 
between the field energy flux and the field momentum densities 



(2i 



W. Pauli refers to the identity as a theorem and makes a comment similar to that 
of C. Lanczos, IPauRI] p. 125]: "This is the theorem of the momentum of the energy 
current, first expressed by Planck^^^ according to which a momentum is associated with each 
energy current. This theorem can be considered as an extended version of the principle of 
the equivalence of mass and energy. Whereas the principle only refers to the total energy, 
the theorem has also something to say on the localization of momentum and energy. "There 
is an intimate relation between the concept of particle and the field concept of the symmetric 
energy-momentum tensor (EnMT) T^^ . In particular, the fundamental Einstein mass- energy 
relation Eg = mgc^ can be interpreted as the symmetry of the energy-momentum tensor, a 
point stressed by C. Lanczos, ILan VPM\, p. 394]: " It was Planck in 1909 who pointed out 
that the field theoretical interpretation of Einstein's principle can only be the symmetry of 
the energy-momentum tensor. If the Ti4 (i = l,2,3j (i.e. the momentum density) and the 
T^i, the energy current, did not agree, then the conservation of mass and energy would follow 
different laws and the principle m = E could not be maintained. Nor could a non-symmetric 
energy-momentum tensor guarantee the law of inertia, according to which the centre of mass 
of an isolated system moves in a straight line with constant velocity. " 

The Noether theorem then implies 10 conservation laws: the energy and momentum 
conservation laws 

dC 

d,T^^ = r, r = -Q^^^ (29) 
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where f is the force density, and the angular momentum conservation law, \Sexl\ Sec. 10.2], 
[BiHItl Sec. II.l, III.4], 

d^M^"'"' = AT^T, where (30) 

]\^tj.u'y jg angular momentum density tensor, and N'^'^ is the torque density tensor. Notice 
that the angular momentum conservation follows readily from the energy and momentum 
conservation law (12^ combined with the symmetry of the EnMT T'^'^. 

Using the interpretation fl271) we can recast the energy-momentum conservation laws fl29|) 
in a more appealing form 

dtp^ = J2 - f'=o-.^' = ^' 2, 3, (31) 

3=1,2,3 

d,u = -J2 d^s^ - f = (32) 

i=l,2,3 

In the case of a closed system when /* = 0, the total conserved quantities are, |Barutl 
(3.76)-(3.77)], IPauRFTl (6), (12), (15b)], Sec. 10.2], 

pu ^ /"T^-da^, = [n'^'^da^, (33) 



where a is any space-like surface, for instance xq = const. is the four-vector of the total 
energy-momentum and J'^'^ = —J'^'^ is the total angular momentum tensor. Importantly, 
for closed systems the conserved total energy-momentum P'^ and angular momentum J'^"' 
transform respectively as 4-vector and 4-te-nsor under Lorentz transformation, and that is 
directly related to the conservations laws, |Moll Section 6.2], |Jackt Section 12.10 A]. But 
for open (not closed) systems generally the total energy-momentum and J^^ angular 
momentum do not transform as respectively 4-vector and ^-tensor, [Moll Section 7.1, 7.2], 
[Jackl Section 12.10 A, 16.4]. 

The formula for the 4-microcurrent density turns into 

ju ^ _Xq_. f ^gu*^* _ ^*Qu^ _ 
2m V / 

Using the same method as in |BF4l Sec. 11.7] we find that for the Lagrangian (fTO|) the general 
conservation laws ( l29l) turn into 

d^Tf-" = /^ (34) 

where 

r = h^F:^ = (/°, f) = Q J ■ Eex, pEe. + i J X Bex) (35) 



with J,p given in (|22|) . Evidently, f is the Lorentz force density for the external field F^^. Us- 
ing the interpretation of the EnMT entries ( l28l) and (!25l) we get the following representations 
for the energy and the momentum densities: 



u = — 
2m 



2 r 



1 



(36) 

{p\p\p') = {dt^pVr + d:r^^) ■ m 
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2.2 Relativistic dynamics of the energy center of a localized charge 

Here we derive equations for the energy center r based on the conservation laws introduced in 
the previous section. Then based on these equations under asymptotic locahzation assump- 
tions we derive the relativistic point mass equations ([T]), ([3]) for the energy center r. We set 
first the total energy, momentum and the force expressions respectively by 

E (t) = J u {t, x) dx, P' {t) = j p' {t, x) dx, F' (t) = j f {t, x) dx. (38) 

The coordinates r* of the energy center r are given then by 

^ x'u (t, x) dx. (39) 



E(t) 



The continuity equation fl23|) when multiplied by (x — r)* readily implies the following ex- 
pression for the current density J: 

dt ((x - r)' p) + pdtr' + V ■ ((x - r)' j) = J\ (40) 

Integrating over the entire space the energy-momentum conservation laws (pTl) and (!32l) we 
get 

^dtE = F°, dtP' = F\ (41) 
For the energy component in fHTl) using ( HUj) we obtain 

F'= [-J- dx = -p-dtr ■ Ee. (t, r) + 6% (42) 
J c c 



'^^ = ^ E / {^oA ((x - r)^ p) + E:,V • ((x - r)^ j) ) dx+ 
+ ^dtr ■ J p (E,, (t, x) - Eex (t, r)) dx. 



(43) 



Integrating the angular momentum conservation law fl30|) for u = and 7 = i over the entire 
space we obtain 

^dt {chP' - r*E) = ctF' - F\' - 6} (44) 



c 

where the particle discrepancy terms 5y are of the form 



5} = J {x'-r') f dx. 



(45) 



The identities dH]) and dB]) imply that 



P^ = \dtr^--5), ^ = 1,2,3. (46) 
c^ c ■' 

Now, combining (H6il with the second equality in ( HTl) . we obtain the relation 

dt(\£dtv--5^=¥, = 1,2,3, (47) 



c^ c 
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with F = defined by (l38l) . (l35l) . The expressions for F can be written in the 

form 

F = fLor(t,r) + (5^, (48) 
where fLor {t, r) is the Lorentz force ([2]) with the remainder term 



j (Eex (t, x) - Ee, (t, r) + c-i^^r* x (Be. (x) - Bex (r)))pdx+ (49) 
+ y"^(^5,((x-r)p) + ^a,((x-r)j')j 



X Bp^dx. 



Equations (H2l) and ( 147|) result in the following system of two ergocenter equations: the spatial 
part 

dt (^^dtv -Sf^= fLor {t, r) + 6f, (50) 

and the time part 

-9tE = iair-fLor(t,r) + 4. (51) 
c c 

Observe now that the particle discrepancy terms <5/,5i?,5^ involve in their integrands the 
factors (x — r), Bex (t, x) — Bex (t, r), Eex {t, x) — Eex (t, r) which vanish at the energy center 
r, and they are small in a small neighborhood of r. In addition to that, the integrands involve 
the factors 3,p which cannot be large outside a small neighborhood of the energy center if 
the solution ip is localized. 

Under the assumption that the charge is localized we neglect the particle discrepancy 
terms <5/, Sp, 5% obtaining the following limit system for r and E: 

dtE = dtv ■ fLor {t, r) , (52) 

dt (^^dtr^ =fLor(t,r). (53) 

Equation ( !53|) evidently has the form of the relativistic version of Newton's law of motion 
with the Lorentz force, namely 

9i(M9ir) = fLor(t,r), (54) 
provided the mass M is given by Einstein's formula 

M = -I- (55) 

Equation f l52|) has the form of the time-component for the relativistic point dynamics, see 
|PauRTl Section 29, 37], |Barutl Section III]. Let us derive now a relation between the 
inertial mass M of moving charge and the rest mass. Using f l52l) we readily obtain 

Mdtr ■ dt (Mdtr) = MdtV ■ qE^^ {t, r) = Mc^dtM, (56) 

which implies the relation 

M2 - {dtvf = Ml (57) 

c^ 
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where Mq is the constant of integration. Consequently, we recover the well-known formula 

M = 7Mo, ^={l-{dtvf lcY^'\ (58) 

The relations and readily imply the accelerated motion equation ([T]). As to the 
equations ^% and they describe asymptotic behavior of the energy and the ergocenter 
of the charge when its wave function if) remains localized in the course of motion. 



Remark 2 We would like to stress again that the rest mass Mq in our treatment is not a 



prescribed quantity, but it is derived in (51) as an integral of motion (or, more precisely, an 
approximate integral of the field equation which becomes precise in an asymptotic limit). As 
any integral of motion, it can take different values for different "trajectories" of the field. This 
is demonstrated by the different values of the rest mass Mq for different rest states constructed 
in the following Section lKM The integral of motion Mq can be related to the mass mo of one 



of resting charges considered in Section \2.4\ by the identity 



Mo = mo (59) 

if the velocity vanishes on a time interval or asymptotically as t ^ — oo or t —)■ oo. // the 
velocity dtV vanishes just at a time instant t^ it is, of course, possible to express the value 
of Mq in terms of E = E {ip) by formulas / f53]) . but the corresponding ip = ifj [to) may 
have no relation to the rest solutions of the field equation with a time independent profile 
. It is also possible that dtr never equals zero, and in fact this is a general case since 
all three components of velocity may vanish simultaneously only in very special situations. 
Hence, there is a possibility of localized regimes where the value of the "rest mass" Mq may 
differ from the rest mass of a free charge. In such regimes the value of the rest mass cannot 
be derived based on the analysis of the uniform motion as in Section \2.4\ This wide variety 
of possibilities makes even more remarkable the fact that the inertial mass is well-defined and 
that the Einstein formula Oj holds even in such general regimes where the standard analysis 



based on the Lorentz invariance of the uniform motion as in Section 2.4 does not apply. In 
a general case where the localization is not assumed, the functional 

^0 = ^ (l - ^ (dtrf^ (60) 

extends formula (5^ to general fields and produces a "generalized rest mass" Mq defined in 
terms of the energy and the ergocenter as follows: 

Mo2= hm i- r Ml{t)dt. (61) 



T^oo 2T 



T 



The above formula obviously defines the value of the rest mass for a more general class of 
field trajectories, and according to ( f53]) and (51) produces the value of the integral of motion 
in the case of asymptotic localization. 

Remark 3 An accelerating particle in an external field is not a closed system, and there are 
principal differences between closed and non-closed systems. In particular, the total momen- 
tum and the energy of a closed system are preserved. For closed systems the particle equation 
and the momentum kinematic representation can be derived from the field theory with the use 
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of the angular momentum conservation, lMoI\ Sec. 7.1, 72], ISexH Sec. 10.2], lLanVPM\ Sec. 
23]. For non-closed systems the center of energy (also known as center of mass or centroid) 
and the total energy-momentum are frame dependent and hence are not ^-vectors, \Mol\ Sec. 
7.1, 7.2], ILan VPM\ Sec. 24]. The rules of the transformation for the energy-momentum are 
due to Einstein and Laue, \PauRT\ Sec. 43]. 



2.3 Rest states, their energies and frequencies 



We suppose for a resting charge {ip (t, x)| to be time independent. Such resting states of the 
charge exist in the absence of external fields, ip^^ = 0, Aex = 0. We are particularly interested 
in rest states ijj that vary harmonically in time and consider solutions to the KG equation 
in the form of a standing wave 



(62) 



where ip (x) is central-symmetric. The substitution of fl62l) in the KG equation ([5]) yields the 
following nonlinear eigenvalue problem 



(63) 



Recall now that the solution ip must also satisfy the charge normalization condition (|24|) 
which takes the form 



ax = — , 



mc 



X 



(64) 



The energy defined by ( l36l) . ( I38l) yields for a standing wave ( l62l) the following expression 



X 

2m 



dx. 



(65) 



The problem ([HSD, flM|) has a sequence of solutions with the corresponding sequence of fre- 
quencies oj. Their energies Eqcj are related to the frequency u by the formula 



Xa;(l + e(u;)), 



ac 



X 

mc ' 



(66) 



where ac = |^ is the reduced Compton wavelength of a particle with a mass m and from now 



2-K 



on we assume that x = ^■ 

The above expression for the energy involves the size parameter a. It coincides with 
the size of a free electron in the absence of EM fields. If our electron is placed in a strong 
EM field, its wave function ip significantly changes. In particular, in the hydrogen atom the 
Coulomb field of a proton causes ip to become a perturbation of an eigenstate of the linear 
hydrogen problem, [BF6j . Therefore, its actual size in the hydrogen atom is of order of the 
Bohr radius ob. It is physically reasonable to assume that the Coulomb field of the proton 
causes the distributed charge of the free electron to shrink, therefore a should be larger 
than Ob. An analysis shows that for quantitative agreement with the classical hydrogen 
spectrum the value of a for the electron should be at least IOOob, |BF6] : therefore, in ( 166|) 



a a\ 



\la} ^ 10 ®. This relatively large size of a contrasts sharply with the concept 



of a point charge, but agrees well with physical properties of electron both at atomic and 
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macroscopic scales. Namely, for a larger value of a the nonlinearity becomes smaller, and 
the charge distribution at atomic scales becomes closer to the De Broglie wave, |BF4] . Note 
also that the electrostatic potential generated by the free charge with the size a is very close 
to the exact Coulomb potential at a macroscopic distance R from the center of the charge, 
namely the difference is of order e~^^/"^ and is extremely small if R is greater than a. 
The Gaussian wave function 

(t, X) = e--0*a-3/2^-3/4g-|.|V2a2 ^57) 

with w = Wo is the ground state to the problem (it is referred to as gausson in |Biaj ). This 
state has the minimal energy among all functions satisfying f lM|) . hence it is stable. In f lBB]) 
a':$> ac, and u > uq since loq corresponds to the ground state. Therefore, the relation (E^ 
differs only slightly from the Planck- Einstein energy-frequency relation E = hu. Note that 
in the non-relativistic version of our theory in ( |BF6] ) the relation AE = hAu is an exact 
identity for hydrogenic atoms. 

By a change of variables the original nonlinear eigenvalue problem fl63p can be reduced 
to the following nonlinear eigenvalue problem with a logarithmic nonlinearity with only one 
eigenvalue parameter ^ and a parameter-independent constraint: 



V'^, = G[(^\ij\'^ij,~^^„ j |^i|'da; = l. (68) 



The parameter ^ is related to the parameters in (163|) by the formula 



The eigenvalue problem (l68ll has infinitely many solutions (^„, V'ln)) = 0, 1, 2..., representing 
localized charge distributions. The energy of n > 0, is higher than the energy of the 
Gaussian ground state which corresponds to ^ = = has the lowest possible energy. 
These solutions coincide with critical points of the energy functional under the constraint, 
for mathematical details see |Caz] . [BerLioI] . |BerLioII] . The next two values of ^ for the 
radial rest states are approximately 2.17 and 3.41 according to |Bialj . Putting in (169|) the 
values ^ = ^„ we can find the corresponding values of yielding for ^ the following 
approximate formula 

^-l-f||. (70) 

The difference of the energy of the higher states and the ground state energy is small, it is 
of order ^o^^ = ^^^C-Roo where hcRoo = mc^^ is the Rydberg energy. Since we assume 
that ob/ci ~ 10~^ the difference is comparable with the magnitude of the fine structure in 
the hydrogen atom which is of order mc^a"^. This comparison shows that the cohesive forces 
generated by the nonlinearity are relatively small. 

Note that using the Lorentz invariance of the system one can easily obtain a solution 
which represents the charge-field moving with a constant velocity v simply by applying to 
the rest solution [ilj,ip,0) the Lorentz transformation (see jBF4] . |BF5] and the following 
Section [23D. 

Remark 4 The rest states of higher energies corresponding ^ = with n > are unstable. 
Since the charges are coupled via electromagnetic interactions, there is a possibility of the 
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energy transfer from charges with higher energies to the EM field, making such states very 
improbable in normal circumstances. It is conceivable though that in a system consisting 
of very many strongly interacting charges, such as, for instance, an astronomical object, a 
significant quantity of such states may be present contributing to the total energy and mass 
of such a system. 



2.4 Uniform motion of a charge 

Consider now a free motion of a charge governed by the KG equation ^ where the external 
fields vanish, that is ip^^ = 0, Acx = 0. Since the KG equation is relativistic invariant, the 
solution can be obtained from a rest solution defined by (!62l) . (!63|) by applying Lorentz trans- 
formation as in |BF4j . |BF6] . Hence, the solution to the KG equation ([5]) which represents a 
free particle that moves with velocity v is given by the formula 

^|J (t, x) = ^f,,, (t, x) = e-(^-*-''-'')z^ (x') , (71) 

with ^ (x') satisfying the equation fl63|) . and 

(x') = (x') = a''/% (x'/a) , (72) 



where 7 is the Lorentz factor 



(7-I) V 

— (v-x)v-7vt, k = 7cu— , (73) 



l={l-P')-''\f3=U. (74) 



All characteristics of a free charge can be explicitly written. Namely, the charge density p 
defined by the relation (122|) and the total charge E equal respectively 



p = iq 



2 



^ (x') , p= I piji^) (bL = q, (75) 



E = 7mc2 (1 + (w)) , (76) 

where 0(0;) is given by fl5^ . The current densityJ, the total momentum P and the total 
current J for the free charge equal respectively 

^ I /|2 k1 ^ ^ ^ 

m ip m & 



J = ^nim-f IV'I' = 7/if:a;^v (x') , (77) 



P = 7mv (1 + (tu)) , J = y J (x) dx = gv. (78) 

The 4- vector (E, P) is a relativistic energy- momentum 4- vector with the Lorentz invariant 
E^ — c^P^ = (1 + 0)^m^c'^. Hence, based on the commonly used argument, |MoH Sec. 3.3], 



|PauRTt Sec. 37], it is natural to define the rest mass of the charge in terms of the mass 
parameter m by the formula 

mo = m(l + 0(w)). (79) 

A direct comparison shows that the above definition based on the Lorentz invariance of 
uniformly moving free charge is fully consistent with the definition of the inertial mass which 
was derived from the analysis of accelerated motion of localized charges in external EM field 
in the previous subsection, see Remark |2] for a more detailed discussion. 
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2.5 The spectroscopic and inertial masses 

Our theory produces a description of the hydrogen atom, |BF6] . |BF7] . yielding in the non- 
relativistic case an asymptotic formula En = —hcRoo/n^ for the hydrogen energy levels with 
the factor i?oo given by the formula 

In the relativistic case a more complex formula can be derived which involves an equivalent of 
the Sommerfeld fine structure with the same factor hcR^o- The constant i?oo in f lHOj) coincides 
with the expression for the Rydberg constant if m = mg is the electron mass and q equals 
the electron charge. It seems, therefore, natural to refer to the mass parameter m as the 
spectroscopic mass. We make a distinction between the spectroscopic mass and the inertial 
mass since in our theory the mass parameter m of a charge is somewhat smaller than the 
inertial mass mo defined by the formula (!79|) with uj = ujq, namely 

mo = m 1 + --| , ac = . 81 

\ 2 J mc 

The difference mo — m depends evidently on the size parameter a. The question stands 
now: is there any experimental evidence which shows that the inertial mass mo and the 
spectroscopic mass m may be different as in our theory? 

The quantum mechanics and the quantum electrodynamics allow to interpret the spec- 
troscopic data and extract from it the mass of electron known as the recommended value of 
the electron mass m = me = Aj- (e) . The recommended value of the electron mass A,- (e) 
when expressed in units u, can be found in |Mohr6t Table XLIX, p. 710]: 

A, (e) = 5.4857990943 (23) x 10"^ [4.2 x 10~^°] (recommended value), (82) 

where, we remind the common convention, 5.4857990943 represents the mean value of the 
experimental data, 0.0000000023 represents its standard uncertainty (deviation) and 4.2 x 
IQ-w represents its relative uncertainty. The value A^ (e) is found based on very extensive 
spectroscopic data. 

But there is also another class of measurements, namely the Penning trap measurements, 
which can be considered as the most direct measurement of the electron mass as the inertial 
one, and it gives the following mass value, |Farj : 

A, (e) = 5.485799111 (12) x 10"^ [2.1 x 10"^] (Penning trap). (83) 

Observe now that the value of the electron inertial mass coming from the Penning trap 
measurement ( !83l) is larger than the spectroscopic mass ( l82l) . and the difference is statistically 
significant. Indeed, using standard statistical analysis we obtain for the difference 6Ay. (e) = 
mo — m the following 

mo-m = 5A,{e) = 0A7{12) X lO-^\ [3.1x10"^]. (84) 

If we take the recommended value A^ (e) from (182|) as our "golden standard" then the mean 
value of the Penning trap measurements from ( |83l) is more than 7 standard deviations above 
the mean recommended value A,- (e) = 5.4857990943, and the same is true for approximately 
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one half of all Penning trap measurements which are greater than the mean value. This simple 
statistical calculation makes our point, namely that the inertial mass mo as represented by 
Penning trap measurements (E^j is larger than the recommend value m from [8^) . 

We have not succeeded in finding in the existing literature an explanation of the above 
mass difference, but in our theory such a difference including its positive sign can be easily 
explained by the finite value of the electron size parameter a. The resulting relation between 
the two masses is given by the formula f lHT]) which can be recast as 

-| = 2^ . 85 

From relations ( 184|) and ( 182|) we obtain approximate inequalities 

0.87 X 10^9 ^ rno-rn ^ ^ ^ ^^^^^ 
m 

If we assign the effect to the non-zero value of a, using the formula (l85ll we obtain 

0.97 X 10^ ^ — ^ 2.4 X 10^ (87) 

These inequalities are consistent with our prior assessments for a to be at least of order 
10^a~^ac. The mentioned assessments are made in |BF5j . |BF6j based on the analysis of the 
frequency spectrum of the hydrogen atom, see also a discussion after relation fl66|) . Note that 
in the derivation of flHTl) we used an asymptotic formula for the spectrum from our theory 
of the hydrogen atom. Consequently, the terms which were neglected may result in a much 
larger interval for a/ac, and, hence, the relations ( 187|) are more for illustration purposes 
rather than for an actual estimate of the ratio a/ac. 



2.6 Charge localization assumptions 

The ergocenter obeys the relativistic version of Newton's equation if the particle discrepancy 
terms <5/, dp, 5^ in the ergocenter equations (15T1) . ( I50ll can be neglected. These terms would 
vanish exactly if the charge and corresponding currents are localized exactly at the center r, 
or if r is a center of symmetry and the external EM potentials are constant as in the case of 
a uniformly moving charge. In the general case we may only expect that these terms vanish 
asymptotically in a certain limit, namely, 

8f ^0, 8f-^ 0, 4 ^ 0- (88) 

There are two kinds of quantities which enter the discrepancy terms. One kind involves 
differences Eex (t, x) — Eex (t, r) and B^x {t, x) — Bex {^i r); they vanish if the fields are constant 
and are small if the fields are almost constant. The magnitude of inhomogeneity of the EM 
fields can be described by the typical length Rf at which they vary significantly. Hence, 
to ensure that the fields are almost constant near the charge, we assume that Rf is much 
larger than the charge size a and impose on the external field the following asymptotic local 
homogeneity condition: 

a/Rf 0. (89) 
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Another kind of quantities which enter the discrepancy terms involve factors (x — r) as in 
dp in dUD: 



X Bcxdx 



|i|^9,((x-r)p) + ^9,-((x-r)j')j 
= ^dt j p (x - r) X Bexdx-^ j p(x - r) x StBex^x + j ^ J - <9;Bexdx. 

These quantities vanish for spatially constant Bex, Egx if p and every component of J is 
center-symmetric with respect to r. Similar quantities are present in S^p and dj defined by 
fl45p . For example 



6) = / - (pEcx + ^ J X Bex^ d. 



(90) 



vanish under the same central symmetry assumption. Hence, to satisfy the charge localiza- 
tion condition (IHHj) . it is sufficient to impose two separate requirements: (i) the asymptotic 
condition (189!) and (ii) asymptotic central symmetry of p and components of J. 

To clarify the meaning of the above localization conditions, let us look at the simplest case 
where p and J are derived for the uniform motion of a free particle described by the solution 
ip (t, x) = ipfj-ee ('^5 ^) considered in Section 12.41 The solution has the following properties: 

(i) the energy density u is center-symmetric with respect to r (t) = vt, hence the ergo- 
center coincides with r (t); 

(ii) the charge density p is given by fl75l) . and according to fl72|) it converges to q6 (x — r) as 
a — > where 5 (x ) is the Dirac delta-function; 

(iii) the current J is given by fl77|) . its components are center-symmetric and converge to 
the corresponding components of qv6 (x — r). 

Hence, the locahzation assumptions ( 188|) are fulfilled for p and P derived for ^£j,gg {t, x) for 
general fields Eex, Bex which are regular near r when a/i?/ — >■ 0. If the motion is almost 
uniform, namely if the external fields are not too strong, and the solution ip {t, x) of ([5]) is close 
to ^/'fi-ee (t, x), we may expect that such a solution also satisfies the localization conditions, 
we present an example in Section [31 

Now let us briefly discuss condition (l89l) . Suppose that the corresponding forces, ac- 
cording to fl53l) . are of the order /Lor where /Lor is a typical magnitude of the Lorentz force 
fLor (t, r). The spatial scale Rf at which the forces associated with the electromagnetic fields 
Eex, Bex vary by the same order of magnitude as /Lor can be defined as follows: 

^ ^ max ( - iVBexl + iVEexl ) (91) 



Rf /l„, k-r|<ea VC 

with 6^1. In view of this definition, condition ( |89ll ensures that the variability of EM fields 
causes a vanishing perturbation to the Lorentz force in (l53l) . 

The estimates of smallness of the discrepancy terms can be made in certain asymptotic 
regimes, but they are laborious, and we would like to make some guiding comments. The 
KG equation fl93l) evidently involves five parameters, namely, c,h,m, q,a, and the external 
fields also can depend on parameters. The limits ( 188|) can be considered simultaneously with 
certain combinations of the parameters tending to their limits, and the limits in f l55]) have 
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to be taken together with these parameter hmits. Note that all the quantities which enter 
equations (!50|) (15T]) are the integrals of certain densities over the entire space. The integrals 
over the entire space obviously are the limits of the integrals over the domains |x — r| < 9a 
with ^ — )■ oo, and the value of 9 in the asymptotic regimes can be related with the values of 
the other parameters mentioned. An example of such asymptotic situation with nontrivial 
relations between the parameters is given in the following section. 

3 Relativistic accelerated motion of a wave- corpuscle 

Up to now, we were considering relativistic features of an accelerating charge assuming its 
localization. The localization assumptions, though natural, are rather technical when it comes 
to rigorous treatment. With that in mind, we study in this section a particular case where 
the dynamical problem is non-trivial and produces a wide variety of accelerated relativistic 
motions of the charge which are simple enough for a detailed analysis and the verification 
of the localization assumptions. The analysis is still rather involved, and that comes at 
no surprise since for general external EM fields the KG equation ([5]) has no closed form 
solutions. We succeeded though in finding a large family of non-trivial regimes for which 
we obtain almost explicit representations of solutions allowing for a detailed study of the 
relativistic features of the charge accelerating in an external EM field. 

3.1 Rectilinear charge motion 

In the previous section we studied dynamics of localized accelerating waves in general EM 
fields. The relativistic point dynamics was derived under the assumption that the localiza- 
tion conditions f l55]) hold. Here we present an example of an accelerating charge when its 
localization can be maintained in the strongest possible form. Namely, we consider a regime 
when the accelerating charge exactly preserves its Gaussian shape up to the Lorentz contrac- 
tion; we call such a solution of the KG equation a wave-corpuscle. The charge wave function 
is similar to ■i/'fj.gg (t, x) in (17T]) and consequently it is localized about r(t), but its velocity 
V = dtv is not constant and the charge has a non-zero acceleration. The shape of the wave in 
an accelerated regime is exactly the same as for a free particle, and it is only the phase factor 
that is affected by the acceleration caused by the external force. Such an accelerated motion 
with a fixed shape is possible only for a properly chosen potential ip^^. We consider here the 
simplest but still non-trivial case where the charge moves and accelerates in the direction of 
the axis with the potential y^gx being a function of only the variable x^ and the time t, 
and there is no external magnetic field, that is Agx = 0. 

When the external potential v^ex depends only on t, X3, the equation ([5]) in the three 
dimensional space with a logarithmic nonlinearity ( IT2|) can be exactly reduced to a problem 
in one dimensional space by the following substitution 

^ = n-'/^a-' exp (^-^ (x? + x^)^ ^j,^ {t, X3) , (92) 

We obtain a reduced ID KG equation for ip = ip^j^ (tyX^) with one spatial variable: 

- \dA^l^ + dli^ - Gil ^ - = 0. (93) 
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where 

^0 = dt=dt + -y^ex, fe^ = ^cx 2:3) , 

and the ID logarithmic nonhnearity has the form 

G'ai im') = -a'' [In (7rV2 + i] _ a-'lna. (94) 
This equation has a Gaussian as a rest solution, 

From now on we write x instead of X3 for the notational simplicity. 

We consider the center location r (t) to be a given function. We also assume that the 
motion is translational. Namely, the solution [t, x) to fl93|) has the following Gaussian 
form: 

= a-i/V'^«(*'^)V'(a"So^/). y = x-r{t), (96) 
^(t,z) = 7r"i/V'/'7re-^'/', 

where 7 = 7 (t) is the Lorentz factor, 7o = 7 (0), the phase Sq (t, y) has to be determined. 
We look for such a potential (p,,^ that the wave if) defined by fl^^ is an exact solution of the 
KG equation fl93|) . If such a potential is found, the function ip (t, x) describes an accelerating 
charge with the strongest possible localization, namely, a charge with a fixed shape {ipl and 
arbitrarily small size a. 

3.2 Mild acceleration regime 

Below we provide an example of a regime which allows simultaneously localization and accel- 
eration, with a general charge trajectory r (t) subjected to certain regularity conditions. The 
KG equation (!93|) involves five parameters, namely the speed of light c, Planck constant h, 
the mass parameter m, the charge value q and the charge size a. A combination of the param- 
eters which is important in our analysis is the reduced Compton wavelength = Now 
we describe relations between the parameters which allow for localized accelerating charges. 

We begin with introducing a class of admissible trajectories r{t). We assume that 
the charge does not undergo violent accelerations and suppose a regular dependence of the 
normalized velocity (5 on the dimensionless time r = ct/a, namely 

max {\drP\ + \dlP\ + \dlP\) < e, where r = -t. (97) 

We also assume that the velocity v = dtr, as well as its variation, is smaller than the speed 
of light c, namely 

1 1 

-max It; (t)| < ei, - max If (t) - t; (0)1 < ei ei < 1. (98) 

c t c t 

Obviously, the variation of f3 and of the Lorentz factor 7 remains bounded for all r, that is 

|7o-7(T)|<Cei, |/3o-/3(t)| <ei, (99) 
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where 70 = 7 (0) , /Jq = /3(0). Observe that the above assumptions are consistent with 
significant changes of velocity v in asymptotic regimes where c — )■ 00, v — 00, f/c — )■ const. 
Note that condition (!99|) does not follow from (!97|) . It can can be seen from the following 
example of oscillatory motion clarifying the above conditions. 
Let V it) be defined by 

V (t) = c/3o + vi sin^ (^^^) = ^/^o + Vi sin^ (r^r) for t > 0, (100) 
(t) = c/3o for t < 0, 

where vi is the amplitude of the variable part of velocity, 1] is a dimensionless frequency, 
/9q < 1 is the initial normalized velocity. We readily obtain the following estimate of the 
normalized velocity variation: 

^\v{t)-v{0)\ < hi, 
and to satisfy (!98|) and (p7|) we set 

-vi = ei, -vi (r] + 7]"^ + T]''^) = e. (101) 
c c 

Since the parameter rj can depend on ei, the boundedness of ei and e do not follow one from 
another. The above conditions imply a bound on the acceleration, 

\dtv\ = r]Vica~^ |sin (2?7ct/a)| < ec^/a, 

which involves a large factor c^/a; therefore, the acceleration can be large even if the param- 
eter e is small. 

Let us consider now the following important dynamical problem: determine the external 
electric field which generates a purely translational motion described by ( 1961) . The size 
parameter a sets up a natural microscopic length scale. We assume that the velocity v (t) 
satisfies conditions (1971) and (p8|) , and the parameters of the problem satisfy the following two 
restrictions. First, we assume that the Compton wavelength defined in ( 166|1 is much smaller 
than the size parameter a, namely 

C = - « 1. (102) 

a 

The second condition relates ( to the parameters e and ei from (!97|) . ( 198|) : 

^C' < 1- (103) 

For a given trajectory r (t) we find then a potential (p^^ such that the Gaussian wave function 
with the center r (t) is an exact solution of fl93|) in the strip 'E (6) = {\x — r {t)\ < 6a} around 
the trajectory, and 6 grows to infinity as C — ^ 0. The solution ijj is similar to the free solution 
but it allows for an acceleration. We will show that the trajectory r (t) relates to the potential 
according to the relativistic version of Newton's law 

dt (m^v) + qd^^ip^, (t, r) = 0. (104) 

Here ip^,. (x) is the leading part of the potential ip^^ which causes the acceleration of the 
charge, we call it "accelerating" potential, it does not depend on the small parameter (. The 
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remaining part of the external potential (p^^ in the KG equation (193|) is a small "balancing" 
potential 

^2 (t, X, = </?ex (t, X; - f^c (t, X) , 

which allows the charge to exactly preserve its form as it accelerates. The balancing potential 
vanishes asymptotically, that is (t, — )■ as C — >■ 0, and the forces it produces also 
become vanishingly small compared with the Lorentz force qdx'Psjc (^) the strip H {9). The 
potential ip^-^ {t, x) converges to its accelerating part v^ac iii the strip S (9) with 9 ^ oo as 
C — 0. In view of fll04p it is natural to treat as a perturbation which does not affect the 
acceleration, but is responsible for the exact preservation of the shape of charge distribution 
during its evolution. In the asymptotic limit C ~^ we have in (I104p 

d.^,Ar{t))^dxP,At,r{t)). (105) 

The summary of what we intend to fulfill is as follows. For any given trajectory r (t) we 
construct a potential <f^^ which makes the Gaussian wave with the center r (t) to be an exact 
solution of the field equation in the widening strip 2(6*), 9 oo. Since the shape of \ip\ 
is preserved, such a function is localized around r{t). We can show that in the considered 
regimes the typical spatial scale R of inhomogeneity of the constructed electric field tends to 
infinity, R/a — )■ oo, and hence fl89l) is fulfilled. 

The implementation of the outlined above approach is provided below and here is yet 
another look at it. We construct the external field which causes the accelerated motion of 
a charge distribution with a fixed Gaussian shape; the distribution exactly satisfies the rela- 
tivistic covariant KG equation in a wide strip about the wave center trajectory. Consequently, 
it provides an example for the relativistic dynamics of a charge wave function with a fixed 
shape. The possibility of a uniform global motion without acceleration is well-known, see 
Section [23] and [BenForj . The fact that relativistic acceleration imposes restrictions on the 
spatial extension of rigid bodies was noted in a different setting in |Eri82j . 

3.3 Equation in a moving frame 

As the first step of the analysis we rewrite the KG equation fl93p in a moving frame around 
r. We take r (t) as the new origin and make the following change of variables 

X3 = r (t) +y, tjj (t, x) = ip' (t, 2/) , v = dtf. (106) 

The ID KG equation (!93l) then takes the form 

+a|V/ - Gl Wi'') i'' - = 0, 

where v (t) is a given function of time. We write the external potential (p^^ which produces 
the motion in the form 

ip,, (t, x) = if,, (t, y) + (t, y, , (108) 
where the accelerating potential is linear in y 
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and {t, y) is a balancing potential^ which, as we will show, is small. The ID logarithmic 
nonlinearity = G'^^ is defined by We assume that the solution ■ip{t,x) has the 

Gaussian form, namely 

ac 

where we define by an explicit expression 

= a-'/% {a-\y) , (110) 
where 7o is the Lorentz factor defined by (!74l) from (p9|) and 

(t, z) = 7^-l/V-^'/^ (111) 

a = a(t) = iln^. (112) 
2 7 

Using (1^ and (EZl), we observe that a is bounded, namely 

a = 0(ei), 9,a = 0(e), d^a = O (e) . (113) 
Substitution of (11091) into fll07p produces 

— 2 1 ^* + *^^^) ~ + ^ 

C \ C ih c 

-ic?,5+iH5 + i^, + ifi-.a,Vvi.' 

n n J 

+ (dy - idyS + iv'^lf'^' - (|v|/f ) v]/' - = 0. (114) 

We show below that and S are of order (ei + e). To eliminate the leading (independent 
of the small parameter () terms in the above equation, we require that the following two 
equations are fulfilled: 

-9,. + ^-.^^7 = -7^.0, (115) 

and 

mdti-fv) + q(p'^^ = 0. (116) 

Obviously, the expressions in flll4p eliminated via equations flllSp and flll6p do not depend 
on (. The equation 0116P determines the essential part ip'^^^y of the accelerating potential; 
the constant part (pQ (t) of the accelerating potential can be prescribed arbitrarily, since we 
always can choose the phase shift s (t) so that f lllSp is fulfilled. Equation (11161) evidently 
coincides with the relativistic law of motion (I104p . 

The above conditions annihilate in the equation (11140 all the terms which do not vanish 
as ( tends to zero; we show this in the following two sections. As a first step, we simplify 
equation (llOTp using equations (IllSp and (11161) . and obtain the following equation: 

+ (dy - idyS + i—--f) ^' - G', f ) ^' - = 0. (117) 

In the following sections we find small and S which are of order [ti + e). 
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3.4 Equations for auxiliary phases 

In this subsection we introduce two auxiliary phases, and then reduce the problem of deter- 
mination of the potential to a simpler first-order partial differential equation for one unknown 
phase. Solving such an equation can be further reduced to integration along characteristics 
and allows a rather detailed mathematical analysis. Here we restrict ourselves to the simplest 
steps in the analysis, but the introduced construction can be used for a much more detailed 
analysis of the relativistic interaction of the EM field with a rigid charge. 
It is convenient to introduce rescaled dimensionless variables z, r: 

r = -t, z=—y = -y, ^ = a^/^^', (118) 
a ac a 

1 C c c 

dy = -d^ = — d^, dt = -dr = C — dr. 
a ac a ac 

We introduce auxiliary phases Z and $, 

Z = Cd^S, (119) 

$ = -CdrS + C(3d,S + (120) 

cn 

they will be our new unknown variables. These auxiliary phases determine the balancing 
potential which we intend to make vanishingly small. Obviously, if we find Z and we 
can find S by the integration in z and we set S* = at z = 0. After that can be found 
from fll20p . Consequently, to find a small we need to find small Z and $. 
Equation (JTTTj) takes the following form: 

-(C9. + i$-i7-/3C'9.)'* (121) 
+ {(d, -iZ + i/37)^ ^ - CG[ i^*^) $ - ^ = 0. 

We look for a solution of (11211) in the strip S in time-space: 

E = {{t,z) : -00 <t < 00, \z\<e} (122) 

where ^ is a large number. Note that = Tc~^^^e~"'o^'^^'^ is smaller than vr'-'^/^e"^^/^ outside 
S and is extremely small for large 6, and later we are going to make 6 arbitrarily large. We 
expand (I12ip with respect to $, Z yielding 

- i$ {(dr - i7 - f3(d,) ^ - i (Co*, - i7 - f3Cd,) ($^) + (123) 
+iZ {Cdz + i/37) * + i {Cdz + i/37) {Z^) - = 0, 

where we denote by Q\E' the term which does not involve $ and Z explicitly: 

g = ^ (- m - i7 - Kdzf ^ + m + '^hf ^) - c'c; {^>*m) - 1. (124) 

Since f lll2p holds, the imaginary part of Q is zero, 

Im g = 2Cidr(T + Cdrl = 0. (125) 
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Hence Q is real and we have 

Q = ReQ = {-^dl^ + d^fi^d,^ + |/39,9,vl/ + ^^^dl^ - G[ {^')^ (126) 
By (El]) with ^2 = e^'"^?, we obtain that 

Obviously, does not depend on a, and has the form( IllOp . but with cr = in the 
definition of tpi in (lllip . Such a i/ji satisfies the equation 

dlk - Gil k = 0, (127) 

and from (11261) we obtain that 

2 2 2 2 /' \ 

Q = + + ^/39.5.vp + - 4) + 2CV- (128) 

Now we rewrite the complex equation (I123p as a system of two real equations. The real part 
of (I123P divided by \1/ yields the following quadratic equation 

g - 27$ + $2 - 2/37Z - = 0. (129) 

The small solution Z to this equation is given by the formula 

Z = (r, $) = -/37 + - 27$ + PV + Q) li • (130) 

The imaginary part of (I123p divided by yields 

- 2$ {dr - Pd,) In ^ - {dr - PO,) ^ + d,Z + 2Zd, In ^ = 0, (131) 

the coefficients are expressed in terms of the given \E' defined by (IllOp . (Illip . and can be 
written explicitly: 

/ I X -1/2 

djn^ = -^lz, dr\n^ = dra, ^,= h--v^O)j . (132) 

To determine a small solution $ of (113 ip . (I130p in the strip S we impose the condition 

$ = Oif;z = 0, -cx)<r<cx). (133) 

The exact solution $ of the equation (I13ip . where Z = (r, $) satisfies (I130p . is a solution 
of the following quasilinear first-order equation 

dr^ - - 0$ (r, $) 9,$ = -2$C {dr - Pd,) In ^ + 209, In ^ + 0,, (134) 

where 0$ and 0, are partial derivatives of (r, $) (0 depends on z via Q = YieQ) with 
$ satisfying condition (I133p . According to the method of characteristics, we write based on 
(11340 the following equations for the characteristics: 

^ = 1, ^ = -0^(r,$)-/3(r), (135) 

— = -2$ {dr - (3d,) In ^ + 20 (r, $) 9, In ^ + 0„ (136) 

(JjO 



24 



with the initial data derived from f ll33p on the hne z = 0: 



Ts=0 = To, Z,=o = 0, $^=0 = 0. (137) 

Note that \E' = n~^^'^e~^o^^/'^+^ is a given function of t,z, and /3,7 are given functions of r. 
By the method of characteristics, the solution of (11341) . (11331) can be found by integration 
along the integral curves. If the right-hand side in (I136P vanishes at $ = 0, the exact solution 
$ would be zero along the curve for every tq; hence it would be zero in H. Therefore the 
magnitude of $ on the curve in the strip S is determined by the magnitude of the remaining 
term, 

Fi = 2e(r,0)9^1n^ + e^(r,0), (138) 

in the right-hand side. The system (I135p - (ll37p is still too complex to hope for a closed form 
solution, but here we want only to show that this solution is small; this implies, in turn, that 
the potential is small as well. To this end, we show below that (I135p - (ll36p is a small 
perturbation of a simpler system which corresponds to the principal part of (11340 for small 

C- 



3.5 Small auxiliary phases and small balancing potential 

Now we estimate the leading part of $ for small (. We obtain from (I128P the estimate 

\Q\ + \dM + \drQ\ = O {eC + e,C) ■ (139) 

Obviously, ii Q = then 9 (r, 0) = 0. Since in the formula (11300 /3^7^ > 0, this formula 
determines for small ( a smooth function B of small $ (r, z) in the strip H (6) , and one can 
see from (11280 that the expression for Q is a quadratic polynomial in z with the common 
factor ('^ and bounded coefficients. Therefore O (r, $) is a regular function of small $ in the 
strip S (9) as long as 

C9^ + C^«l. (140) 

Hence we can take 

e^oo as C ^ 0, (141) 

and, therefore, the width 6 of the strip where we have the exact solution of (11340 . (I130p is 
arbitrary large. If $ = O (C e -|- ei) (certainly this assumption is true for small 1^1 thanks 
to (I133p . the argument below shows that we can extend this estimate to the entire strip S (9)) 
we have 

Z = e (r, $) = -1$ + O (C'e + eiC') • (142) 
This implies the following estimate of Fi in (I138p 

F, = {('e + . 

If we replace G by its principal term from (I142p . we obtain from (I135p - (ll36p a simpler system 
for the resulting approximation $: 

dr _^ dz ^ 1 /. 
ds ds p 

d^ o 2 o 

— = -2<I)(9, -/39,)ln^- -<l>9,ln^, (144) 
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A general solution for this system can be written explicitly. Namely, the characteristic curves 
are given by 

r = ro + s, z = z{t,To)= [ ^ (l - rfr, (145) 
and the value of the phase $ is given by the formula 

^(-'-) = ^K0)^^. (146) 

The integral curves r, z (r, tq) obviously cover the entire (r, z)-plane. Since the original prob- 
lem ( I135p - fll36p is a perturbation of the above system of order ^^e + C^ei, the nice properties 
of (11431) . (I144p imply that the characteristic equations (I135p - (ll36p and their integral curves 
also have nice properties, and $ and $ are close one to another. According to (I137P we set 
^ (''"o; 0) = 0, hence $ (r, z) = 0, and the exact solution $ of (I134p . (I133p in the strip S is of 
order ('^e + C^ei for small (. 

Let us verify now that the balancing potential is small. The solution of (11340 . (11330 . 
obtained by integrating (I135p - (ll37p . satisfies in E estimates 

|«l>| = O {{C\ + C'ei) |;.|) , + \dl^\ = O {C'e + Ce,) . (147) 

From (fTT9D- (n:20|) we obtain that 

^, = ^0{ee\z\+Ce^\z\), (148) 
q ^ 

and obviously the balancing potential vanishes when C — )■ 0. 

Let us show now that the constructed potential ip^^ = ip^^ + satisfies the asymptotic 
local homogeneity condition fISI?]) . Definition fpTl) takes the form 

_}_ ^ ^^^\x-r\<ea {\dl<pj) ^ max|^|<e(|gfy;,^|) ^^^^^ 
Rf maxt |9a;V5acl a max^ jS^v^^J 

To express the magnitude of the gradient of the potential ip^^ in terms of the parameters 
corresponding to the mild accelerations, we use relations (I116p and (^7^ : 

dz^^ = adyip^^ = -a—dt M = —dr (7/?) , 

resulting in the estimate 

2 

max {dz'p^s.cl = O (e) . (150) 

q 

Using (lligp . (]120p and taking into account the definition of ac, we obtain 

dl^ + drd^Z - pdlZ = -^dl^p,. (151) 

mc^ 

We observe that, since (p^^ is linear, d'^ip^^ = d'^(p2, and using (I147p and (I142p . we conclude 
that in the strip S 

q ^ 
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Consequently, formula (11490 takes the form 




o {e + ^c^) 



(152) 



We obtain that fl89l) holds in the asymptotic regime 



C' -> 



(153) 



which is consistent with the assumptions fllOSp . f ll02p . The above analysis shows also that the 
magnitude of the balancing potential is vanishingly small compared with the magnitude 
of the accelerating potential. 

The summary of the above argument is that we can find a potential ip^^ (t, x) = ip^ + 
such that (i) the function -0' with the Gaussian profile defined by (I109P is an exact solution of 
the KG equation (11070 in a widening strip S {6); (ii) the balancing potential {t, x) vanishes 
asymptotically, and the relative magnitude of the balancing potential (t, x) compared with 
the accelerating potential ip^^ tends to zero as C — ^ 0. The charge distribution ip provides an 
example of a localized solution in its strongest form, namely, an accelerating solution with a 
fixed Gaussian shape. The possibility to preserve localization which we discuss in this section 
concerns microscopic details of the charge evolution, and naturally the ratio ( = ac/a of the 
charge size to the Compton wavelength oc plays an important role. In the context of Section 
12. 6t a is vanishingly small compared with the macroscopic scale Rf as (I152p clearly shows. 
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